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Abstract 

We discuss several aspects of the dielectric response theory application to the 
density functional theory. This field has been an unceasing source of confusion 
during several decades. The most frequent reasons for this confusion are (a) 
uncritical tranfer of the results, especially regarding so-called local field cor- 
rections, obtained in many-body perturbation theory onto density functional 
theory, and (b) mixing up the statements true for the exact density functional 
theory with those applicable to the local density approximation only. In these 
notes we try to draw an appropriate lines between those theories. We also 
discuss a newly introduced (X. Gonze, Ph. Ghosez, and R.W. Godby, Phys. 
Rev. Lett., 74, 4035, 1995) " polarization-Hdensity functional" and show that 
within a given (e.g., local density) approximation to the exchange-correlation 
energy the Gonze et al approach is exactly equivalent to the conventional one. 
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Even now, 30 years after the discovery of the density functional theory (DFT) |T| there 
is still considerable confusion about applicability of DFT to static dielectric response 0. 
On the one hand, there is seemingly no room for questioning the validity of the DFT in 
this respect: The Hohenberg-Kohn theorem, which is the basis of DFT, is a mathematical 
theorem, rigorously proven, that states that the total energy of the ground state of a many- 
electron system is a unique functional of its electronic density, and the density- density 
response function, defined as 

x{v,v')=5^E/5p{v)5p{v'), 

is also a unique functional of the electronic density. Surprisingly, this simple theorem is not 
readily digested by everybody. The reasons for confusion are several: First, one commonly 
confuses the exact DFT with its local, approximate version (LDA). As we shall discuss 
below, mixing up these two notions is much more dangerous when dealing with dielectric 
response than with the total energy itself. Second, to have this theorem satisfied, one has to 
take into account properly exchange-correlation corrections to the dielectric susceptibility — 
and be aware that those corrections are functionally different in DFT, and, say, in the many- 
body perturbation theory or in the Fermi-liquid theory (needless to say that all observable 
quantity are the same in any of these theories). Third, the notorious failure of LDA-DFT 
to produce the correct band gap (well understood by now, thanks to seminal works of 
Levy and Perdew 0] and Sham and Schliiter 0), combined with the conventional wisdom 
that the dielectric gap determines the response, forces ingenuous researchers to question 
the formal applicability of the DFT to response functions calculations. Finally, in the last 
decades accurate LDA-DFT calculations of the dielectric constant of semiconductors have 
been performed which were usually in error by 10-15%. Incredible success of the LDA 
calculations of structural properties (a few per cent) made people imply indirectly that 
similar accuracy is attainable in other ground-state properties calculations. When it turned 
out to be not true, a suspicion rose that something may be principally wrong with the 
approach (in reality, of course, 10-15% accuracy is excellent for such a simple approximation 
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as LDA). In particular, the Hohenberg-Kohn theorem itself has been questioned The 
goal of the current notes is not to produce a new physical result but to help, especially 
newcomers in the field, to avoid confusion and misleading of incorrect claims spread around 
in the literature. 

One of the source of confusion is the fact that exchange-correlation local field (XCLF) 
corrections are formally different in all three theoretical approaches to the dielectric response 
of many-electron systems. The most traditional one is the many-body perturbation theory, 
some times referred to as Green function theory. In this theory one-electron excitations are 
poles of the Green function, G^^{r,r',uj) determined by the Dyson equation 0, 

G-\v, r', uj) = G,\r, r', uj) - S(r, r', lu) 

where S is the self-energy operator, and G^^(r, r', u) is the Green function of free electrons. 
An alternative form for the Dyson equation is 

i/o(r)V'n(r, E^) + J S(r, r', K)^n(r', = K^„(r, K), (1) 

where Hq is the one-electron Hamiltonian, V'nlr, En) is the Green function amplitude and En 
is the corresponding pole. En are in general complex. In the random phase approximation 
the dielectric susceptibility of the system can be written in terms of the bare polarization 
operator as 

xTpAi^,2) = [l-rc;(l-3)n^^(3,4)]-%«^(4,2), 

where we introduced superscript GF for the Green functions theory, and integration over 
the space coordinates is implicitly implied. Here 1, 2, etc are short for ri, r2, etc. Uq is the 
bare polarization operator. Symbolically one can write this equation as 

ttGF 

GF _ ^^0 (n\ 

and we shall use such notation throughout the paper, when unambiguous. Going beyond 
RPA demands replacing 11^^ by renormalized (full irreducible) polarization operator fl"-^^. 
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The XCLF, as introduced by Hubbard in 1958 [0, assumes that 11*^^ can be approximately 
written as 



^^^-u^^/ii-mr), (3) 

where exchange-correlation interaction J(r, r') is negative (as well as Hq^), and thus the 
total susceptibility is enhanced compared to RPA: 

ttGF gf 

^ l-{Vc + mr I-IxWpa 
An important fact is that the difference between 11*^^ and comes first of all from 
exchange processes, namely 





In semiconductors, the exchange interaction remains long-ranged. In other words, when 
this series is approximated by Eq. ^, 



+ ... 



the corresponding effective interaction / (wavy line) diverges in reciprocal space as 
As we shall see below, it is not the case in LDA. 

Formula (^ looks very similar to the expression for the susceptibility in the Fermi-liquid 
theory 

However, while the left-hand sides of Eqs. || and || are the same, Xrpa differs from x'rpai 
as it usually defined in the many-body theory, in the sense that it is calculated with exact 
one-electron excitation spectrum, given by the poles of the full Green function (|1|), while 
x'rpa assumes bare Green functions. For instance, in the theory of the homogeneous elec- 
tron gas x'bj'a is j^st Lindhardt susceptibility. However, in practical calculations for 



semiconductors, for instance in the GW approximation, it is common to call "RPA" sus- 
ceptibility the quantity calculated with fully renormalized one-electron spectrum, that is, 
Xi?pA- Correspondingly, effective interaction /(r, r'), the effective Landau interaction, is the 
analog of Coulomb + Hubbard interaction in the Hubbard approximation. 

The most common way to calculate the dielectric function is related to the density 
functional theory. To refresh readers' memory, let us remind the basic equation of this 
theory: 

HDFT{r)ipn{r) = Ho{r)Lpn{r) + VKsVn{r') = £„V3„(r), 

where the wave functions (p and the (real) eigenvalues e are for a fictitious system of non- 
interacting fermions with the same density as the electronic system in question. The above 
equation is a tool to find this spectrum, and the Kohn-Sham potential Vks is defined as 

^^^^'^ = 6p{r) ■ 

Eint is a unique functional of the total density, and the total energy of the electronic system 
is, by definition, I]e<^(v^i| ~ + -^mt[p(^)]> where /i is the chemical potential. Since the 

dielectric susceptibility can be defined entirely in terms of the total energy and total density, 

v(r r') = -JMpL (6) 

it is also a unique density functional. Unlike the Green functions theory, DFT allows a 
formally exact expression for x, which follows directly from Eq. ^ and looks similar to the 
approximate Eq.^ in the Green functions theory: 

DFT 

y = ^MA (^) 

^ 1 T ^,DFT ' V / 

-L J'XcXrPA 

where now by definition 



4c (r, r 



6p{r)5p{r') 6p{r') ' 
There is principal difference between I^c behavior in LDA and exact DFT: In the for- 
mer, Ixc{r,r') oc 6{r — r'), while in the latter Ixc{r,r') may be of arbitrary long range. A 



good example (which was communicated to the authors by O. Gunnarsson) is exact DFT 
for an insulator: adding one electron to arbitrary large insulator induce uniform shift of 
the exchange-correlation potential for the whole system (so-called "density- derivative dis- 
continuity", Refs. 0,01), thus making Ixc{'r,r') of infinite range. Since the density-derivative 
discontinuity is always non-zero (because of the exchange-correlation part of kinetic energy, 
see Ref. [Q), that means that in the exact DFT Ixd^) diverges at g — > as S{q). Never- 
theless, the magnitude of the divergency may be arbitrary small. In fact, it appears that 
the DFT expression (^, if the spacial dependence is properly dealt with, yields good results 
(with 10-15% accuracy) even in LDA 0. Again, it is instructive to compare the longe- 
distance behavior of I^c in GF theory, in LDA and in the exact DFT: in GF theory I^^ 
diverges, in insulators, at g — ^ as 1/g^; this statement is not necesserily true in exact DFT, 
where the most divergent term is just a (5-function and nothing can be said rigorously about 
the next terms. Finally, in LDA I^c remains constant at g — * OQ 

The spatial dependence is rather important. Proper treatment of the spacial dependence 
includes Umklapp processes, 

l/<q) = [eq+G,q+G']oo (8) 

Ccq+G' = ^GG' + ■ ■ [Sgg' - (Vc+i.c) n]-^ 

where Vc(q + G, q -|- G') = 47re^(5GG'/lQ + and Ixc is the Fourier transform of Ixc{t^, r'), 
and the tensor dot-products are taken in the right-hand side. All terms here which originate 
from G,G' ^ are called "local field corrections" . An elegant approach which is mathe- 
matically equivalent with Eq.^], and is often called "the Sternheimer equation" , avoids using 
polarization operators explicitly, but instead deals directly with the change of the density 
(see, e.g. Ref. 0, Eqs. 12-16: 



^One can hear occasionally statements that I^J^'^ should diverge as l/q"^ at at small g's; we are 

not aware of any proof of this statement which would not indirectly use an unproven parallel with 
jGF 



5p(q+G)=4 (V^nk|e-'('^+^)'-Pc|5V^nk+,q) (9) 

]<.,n£occ 

where Pc is the projector operator on the conduction (unoccupied) bands. These equations 
are formally equivalent to Eqs.^ and thus are exact within the DFT. Furthermore, if one 
uses the same (approximate) exchange-correlation energy functional, for example LDA, both 
sets of equations (P) and (H) should yield the same result. 

Much of confusion was raised by the fact that the Hohenberg-Kohn theorem applies to 
the toto/ electron density. That is, upon applying a long wave (g — > 0) perturbation, the total 
energy depends not only on the periodic part of the density change, Z]g7^o'^p(q+ G), but 
also on the long wave part, Sp{q + 0). Despite this fact, neither Eqs. H nor Eqs. |^ demand 
supercell calculations with periodicity 1/q. The reason is that the original, unperturbed 
system, which is periodic, bears all information (in linear regime) about the properties of 
the system, perturbed by any external field (as long as it weak enough), including long- 
range perturbation. It was claimed recently by Gonze et al that Eqs. ^ are incomplete 
because the do not treat the long-wave part of the density properly. This is of course a fallacy. 
Interestingly, trying to overcome non-existent incorrectness of Eqs. |], Gonze et al derived 
another equation, but failed to realized that it was mathematically equivalent to Eqs. 
This is worth elaborating. 

Gonze et al suggested in Ref. 0] to use a functional of two functions: the periodic electron 
density, Pper = X^Gt^o p(q + G)e*^'^~^*^^'', and macroscopic polarization P(q). The latter is 
uniquely related to the macroscopic density perturbation, p^ac = p(q)e*'^'' ~ z(qr)p(q) for 
g — 0, so P and pp together uniquely define the total density and therefore the total energy. 
Obviously, one has full freedom to chose with which variables to work, periodic density 
plus macroscopic density, or periodic density plus macroscopic polarization. Having chosen 
the second scheme, Gonze et al derived equations where instead of the change of the total 
exchange correlation potential in 6Hdpt they used only periodic part of this potential. Their 
final expression corresponded to substitution of PcSV^°^\ipnk) in Eqs. ^ by 
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G^o '^P'^PG(r) 

In fact, Gonze et al claimed that previous LDA calculations of the dielectric response in semi- 
conductors were incorrect, because they had used Eqs.p instead of Eqs.|TO|. One can easily 
show, however, that the two sets of equations are formally equivalent: From the definition of 
4, it follows that {6^E^J6P^)6P = q/xc(q, q)5p(q), and [6^E^J6P6p{ci+ G)]6p{ci + G) = 
q/2,.c(q, q+ G)5p(q+ G). One can now notice that the part of the Eqs. ^ which depends 
on the macroscopic component of the exchange- correlation potential, SV^"''^, is nothing else 
but 

G 

which is the same as the second term in Eq. |T0[ Not surprisingly, polarization-f-density 
functional of Gonze et al appears to be equivalent to the original Hohenberg-Kohn functional. 
An unfortunate consequence of this fact is that the hope that density-polarization functional 
of Gonze et al can remedy the above-mention deficiency of LDA, namely the local character 
of the interaction /a;c(r,r'), is futile. New, more advanced approximations to the density 
funcional are needed to improve the results. It is possible, although not garanteed, that 
these approximations will be easier to deal with in the density-polarization formulation than 
in the total-density formulation. A promising routes are generalized density approximation, 
much advanced lately, and truly non-local functionals like weighted-density approxiamtion. 
In this regard, a recent study by Dal Corso et al 0, where a sizable improvement over LDA 
was found for the dielectric constant of silicon, provided that calculations are done at the 
same (experimental) lattice parameter. 
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